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We propose a new type of effective densities via the potential distribution theorem. These densities are for 
the sake of enabling the mapping of the free energy of a uniform fluid onto that of a nonuniform fluid. The 
potential distribution theorem gives the work required to insert a test particle into the bath molecules under 
the action of the external (wall) potential. This insertion work W ias can be obtained from Monte Carlo (MC) 
simulation (e.g. from Widom's test particle technique) or from an analytical theory. The pseudo-densities are 
constructed thusly so that when their values are substituted into a uniform-fluid equation of state (e.g. the 
Carnahan-Starling equation for the hard-sphere chemical potentials), the MC nonuniform insertion work is 
reproduced. We characterize the pseudo-density behavior for the hard spheres/hard wall system at moderate 
to high densities (from p* = 0.5745 to 0.9135). We adopt the MC data of Groot et al. for this purpose. The 
pseudo-densities show oscillatory behavior out of phase (opposite) to that of the singlet densities. We also 
construct a new closure-based density functional theory (the star-function based density functional theory) 
that can give accurate description of the MC density profiles and insertion works. A viable theory is established 
for several cases in hard sphere adsorption. 

Key words: potential distribution theorem, hard spheres, adsorption, effective density, structure, integral 
equation, closure 
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1. Introduction 

Inhomogeneous fluids pose, due to their vast varieties, challenges to the liquid-state theories, 
and are of capital importance in many industrial processes. The nonuniformity in the fluid can be 
produced by a one-body potential: such as the surface potential from the interfaces separating two 
or three coexisting phases (solid/liquid/ vapor phase transitions), electrical fields acting on ionic 
species (the electric double layer), or magnetic fields on ferromagnetic liquids. The conventional 
classical density functional theory (DFT) employed to deal with such systems in the last half 
century has been based on the Hohenberg-Kohn [1J formulation of a grand potential, f2, that is 
expressed in terms of an intrinsic Helmholtz free energy functional i^pi}" 1 ] (IHFEF): 

« = F[P$] + I dfp$(f; H){w(r) - Mo } (1.1) 



where p^ is the singlet (nonuniform) density, w(r) the external (1-body) potential, and //n the 
chemical potential of the bulk fluid. At equilibrium, the grand potential is minimized with respect 
to the singlet pw\ and it yields the Euler-Lagrange equation 



■Pno + pw{f)=Q (1.2) 



where j3 = l/(fcT): k is the Boltzmann constant, T is the absolute temperature. The IHFEF can be 
separated into two parts: the ideal part (which is known: F ld [p] = J drplP (f) [In pw \r)A 3 — 1]) 
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and the excess part F cx . Thus ] = ^ ld + F cx . Recently, we have proposed a star-function 

based density functional theory |5J [3] (s-DFT) which expresses _F CX succinctly as 



= 



&rC^\r)5 Pw {r) + ] i f drdr' C { ^ {r,r')5p w {r)8 Pw {r') 



din 



where S w is the star function [3] defined as the "primitive" of the bridge function B u 

71 

5* = / dr — I dxB w (x) 



(1.3) 



(1.4) 



7o 



and C^\r) is the singlet direction correlation function (1 — DCF), C^\r,r') is the pair direct 
correlation function (2 — DCF). Fq x is the bulk fluid free energy. Subscripts and w indicate 
uniform and nonuniform fluid properties, respectively. 8p w = (pw — Po) is the difference between 
the nonuniform density pw and the uniform density po and j w is the indirect correlation. We note 
that if the bridge function B w can be obtained exactly, the star function will consequently be 



exact from equation (1.4 1, thus the excess free energy F GX in equation (1.3) will also be exact. The 



Euler-Lagrange (EL) equation (1.2) can be expressed as 



p w (r) = pexp -w(f) + C«(r) - C { 1] (r) 



(1.5) 



where is the nonuniform singlet direct correlation function (1 



. - - ,., - ... DCF). Alternatively, it is 

possible to show [2] via functional expansions that equation (1.5) can be written as 



Pw (r) = pexp [-w(r) +lw{r) + B w (r)} 



Thus, we have the equality [3 



C^(f)-C^(f)=-y v> (f)+B w (f). 



(1.6) 



(1.7) 



Equations (1.3) and (1.6 1 constitute the basis of the s-DFT theory. 



In this article, we shall explore two essential aspects of nonuniform fluids: (1) the effective 
density that enables the mapping between the uniform-fluid free energies and the nonuniform-fluid 
free energies; (2) the closure relation between the bridge function and other correlation functions. 
The first task is mediated through the potential distribution theorems [51 [B]; the second is based 
on a successful uniform liquid theory [7] developed in recent years. Historically, there have been 
interchanges and cross-fertilization of uniform and nonuniform liquid theories [H [9] . In this paper 
we test the new formulations on hard spheres adsorbed on hard planar walls. 

In weighted-density density functional approaches (WDA) |10l ITT] the "mapping" of the free 
energy is done through some weighted density: p. The weighted density is obtained from the true 
density P w^ via a convolution integral with a weighting function, to, as the kernel: i.e. p(f) = 
J df 'pw (f')u(r, r'). We shall, however, develop an effective density without going through this 
weighting procedure. This is made possible by the potential distribution theorem (PDT), since 
PDT naturally lends itself to yield the work of insertion (the free energy of an inhomogeneous 
system) . 

The PDT is a fundamental theory of statistical mechanics (5| IB"} [TTHT5] . It is concerned with the 
work required to insert test particles in an equilibrium ensemble representing a uniform fluid or a 
nonuniform fluid. In an earlier paper |14| we have generalized the PDT from the chemical potential, 
to the cavity functions, and to higher-order correlation functions. The PDT has previously been 
extended to nonuniform systems jSJIHI- However, here we distinguish two types of test particles. 
The first type is for a test particle subject to the wall force as well as to the forces of the bath 
molecules. This is the commonly studied type. The second type is for a test particle free of the 
wall force. Its interaction with the wall is intentionally removed (similar to the cavity function 
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in homogeneous fluids). The PDT of the second kind is directly related to the Euler-Lagrange 
equation [3J. The PDT establishes a connection among (i) the insertion work W- lns , (ii) the bulk 
fluid chemical potential /Urj, and (iii) the singlet distribution function pw . To illustrate the PDT, 
we carry out calculations on the system of hard spheres on a hard wall (HS /HW) , for which MC 
data are available (e.g. from Groot et al. |15|). 

The mapping between the uniform and the nonuniform systems can be constructed as follows: 
given an equation of state for the uniform fluid (which is utilized to calculate the uniform chemical 
potential /3/xo), we solve for the hypothetical densities /9 pse udo(z) that will reproduce the insertion 
works W inB (z) obtained from the MC, i.e. /3 Ho(p P seudo( z )) = fiW- ms (z). This density profile is an 
artificial construct and is called the pseudo-density p pse udo(z). Its sole role is to reproduce the 
nonuniform free energy via the bulk equation of state. For instance, in the HS/HW system, the 
bulk fluid equation is taken to be the Carnahan-Starling (CS) [IB] equation which is known to be 
highly accurate. It yields the uniform hard-sphere chemical potential /3[1q. If we know the insertion 
work (3W- ms {z) (from MC) in HS/HW, a series of values of p vsen d (z) can be generated so that 
the equality (i$ s = (3Wi ns (z) holds. The p pse udo(z) obtained thusly will enable the mapping of 
the free energies ( f3p^ s ) of the uniform fluid to those (fiW- ms ) of the nonuniform fluid. In this 
respect, the pseudo-densities perform precisely the same function as the weighted densities p(z) 
in the WDA. The difference is that our approach is based on the PDT, an exact theory, not on 
approximate theories (such as the Percus-Yevick equation |17|). 

Should the bridge function B w be known, one can obtain the fluid structure at the wall via 



the Euler-Lagrange equation (1.6) [2, 3J. B w as a correlation function is well defined in terms of 
infinite series of cluster diagrams |l8| or functional expansions [19]. However, direct evaluation from 
its definition is numerically intractable. But it can be inverted from machine data (as in reverse 



engineering [50]) if we consider equation (1.6 1 as the defining equation |19J . The inverted B w can 
be given as a set of numerical data [5], say, in the form of tabulated entries; or it can be expressed 
as a function B w = f(r) for some analytical function / (e.g. a polynomial). Another approach, and 
a more appealing one is to formulate a liquid theory (a la mode of the Percus-Yevick equation [IT] ) 
that relates B w to other correlation functions B w = f("f) , 7 being a well-defined quantity in liquid 
state theory. This in uniform liquid theory is called the closure relation. This relation may or may 
not exist, because B w is in fact a functional [3JHT] B w = 3>[pw ,w] of the singlet density p$ and/or 
the external potential w ($[.] being the functional sign). Thus, it may not be simply related to any 
one correlation function. For uniform liquids, the determination of the function B = f(j), 7 being 
some correlation, is already a major task (examples such as the Rogers- Young [22\ (RY) closure 
or the Martynov-Sarkisov [23 (MS) closure). We have proposed a theory upon resummation of 
the functional expansions of B w in [3J. This resummation seems to work well for a number of pair 
potentials. For nonuniform fluids, the PY and hypernetted-chain [21] (HNC) closures have been 
tried and shown not to work well in the past [25] . And this prompted the proposition that closure- 
based approaches might not be suitable for nonuniform systems. In this work, we shall propose 
a zero-separation type (ZSEP) closure derived from the study of uniform liquid theory [7 . This 
closure will be tested here on selected cases of the HS/HW systems. Once the singlet density is 



obtained from equation (1.6) given B w , one can use the Euler-Lagrange equation (1.5) to obtain 



the 1-DCF Cw ■ The 1-DCF is simply related to the insertion work as 



pW- w {f) = -C&(f). (1.8) 



Thus, we can access the insertion work through (1.5 1, using the MC data [15] on pffl as input. 
The thusly obtained f3Wi ns is considered as the MC-derived insertion work. In addition, we shall 
construct a theoretical method for obtaining the insertion work. 

Section 2 gives the potential distribution theorems for nonuniform fluids with general interaction 
forces. Section 3 shows the calculations of the effective densities based on the MC-derived insertion 
works for the hard spheres/hard wall system. The mapping of the nonuniform free energies is 
mediated through an accurate equation of state by Carnahan and Starling. We characterize the 
behavior of p pseu do and compare with the commonly used weighted densities found in literature. In 
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section 4 we test the new closure theory (ZSEP) on the same systems. We compare the theoretical 
structures and free energies with the MC data. In section 5 we draw the conclusions. 

2. Potential distribution theorems for nonuniform systems 

In this section, we shall introduce the potential distribution theorems for two types of test 
particles. One type of test particle is the well-known one, i.e. it is subject to interactions with the 
bath molecules as well as to forces of the wall. The second type is the one that interacts only with 
bath molecules and not with the wall. The latter result will be the PDT that we shall use in this 
article. We start with the averages in a canonical ensemble. (This can also be done in the grand 
canonical ensemble |5J|B]). The N-body system consists of TV fluid molecules and they border on 
one side at a solid surface. The Hamiltonian Hm is written as 

H N (p N , r") = K N (p N ) + U N (r N ) + W N (r N ) , (2.1) 

where r N = (fi, r*2, . . . , r/v) is a shorthand for the TV-vector of the positions of TV particles, and p^ 
is the TV-momenta vector, iV/v is the kinetic energy, while 

N-l N 

U N (r N )= E u(23 (^') ( 2 - 2 ) 

l=i<j j=2 

is the total potential energy between fluid particles (assumed to be pairwise additive), and 

N 

W N (f N ) = ^w(r k ) (2.3) 

k=l 

is the sum of one-body energies arising from the external (wall) potential w. This is the potential 
energy that drives the inhomogeneities in the system. 

Next, we consider an (N-l)-body system: the sum of pair energies becomes 

N-2 N-l 

U N _ 1 (r*- 1 )= £ 5> (S, Wi) (2-4) 

1=»<2 J=2 

and the sum of one-body energies 

N-l 

W N - 1 (f N - 1 )=J2 w (rk)- (2-5) 
fc=i 

A test particle is introduced to the (N-l)-body system as particle with label "N" located at 
the distance rjy. This test particle can interact with the other TV — 1 bath molecules through the 
pair potential u^(ri, j-jv), i < TV, as well as with the wall through the one-body potential w(f N ). 
This particle is designated the type-1 test particle. The excess potential energy 'J'™ (r/v) due to the 
presence of type-1 test particle is 

N-l 

* w (r N ) = U N (r N ) - Ujf-!^- 1 ) + W N (r N ) - W^-i^" 1 ) = £ u^(r u f N ) + w(r N ). (2.6) 

i=i 

On the other hand, the type-2 test particle, by choice, does not interact with the wall (w(rjv) = 
0). The excess potential energy ^-^(rjv) is (subscript "-w" stands for "without the wall interac- 
tion") 

N-l 

*- w (r N )=J2 u(2) (ri,?N)- (2.7) 
»=i 
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The potential distribution theorem for the first-type test particle (interacting with w) is obtained 
from the (N-l)-body canonical ensemble average of ^(nv) 



ln(exp [-/3V w ]) N -i ;w = -/3p, w + In [p£ 1 )(f JV )A 3 



(2. 



since 



(exp [-py w ]) N -i y 



{N-l 



dr N - x exp [-PU N -i - PW N - X ] 

'N-l 



x exp < — j3 



uW(r m )+w{r N ) 



i=l 



(2.9) 



Furthermore, from the definition of the singlet density pw\f]^) 

N 



dr^exp [-/3U N -fJW, 



N 



(2-10) 



Equation (2.8 1 has been obtained earlier Note that Qn is the N-body configurational 

integral; A is the de Broglie wavelength; and 



/^ = m%^+ln(pA 3 ). 



IN 



(2.H) 



The potential distribution theorem for the type-2 test particle (without w) is obtained instead 
from the (N-l)-body ensemble average of the wall-less excess potential 9— w (fff) equation (|2.7[), 



i.e. 



ln(exp [-^- w ]) N -i. w = -Pn w + Pw{r N ) + In p$(r N )A 3 

'N-l 

^ u {2 \r lN ) + w(r N 



(2.12) 



(exp [-^_ W ]) N _!. 



Qn- 
x exp < — /3 



'N-l 



pw. 



AT— 1 1 



(2.13) 



The insertion work W- ms required to insert a type-2 test particle with the excess potential 
^-w{?n) is thus (noting that j3p w = /3^ig x + ln^rjA 3 )) 



(3W h 



- ln(exp [-P^- w ])n-1;w = Pfi>w ~ fiw(r N ) - In p^ 1) (r A r)A 3 
/3^ x - fJw{r N ) - In \p${fs)/po 



(2.14) 



The quantity Wins can be directly simulated via the MC method using the Widom particle- 
insertion technique HH], or from the EL equation (1.5|. Wi ns can also be considered as the 
intrinsic work of insertion. 



3. Effective density and work of insertion 

For the hard sphere system, we adopt the MC data by Groot et al. [TS] as the basis of calculation. 
The number densities chosen for the hard spheres are p* = 0.5745, 0.715, 0.758, 0.813, and 0.9135. 
The hard spheres are adsorbed on a planar hard wall stretching over the x-y plane. Thus, the 
fluid-fluid pair potential is 

U( 2 )(r)=oo, r ^ a, 

W 2 )(r)=0, r>a. 1 J 
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Figure 1. The effective (pseudo) density p pBeu do(z) (symbol *; from PDT) and the insertion 
work j3W- lns (symbol A; from MC) at the bulk hard-sphere density p* = 0.5745 for the HS/HW 
system. Symbol ({>) = the singlet density profile pw\z) from MC data of Groot et al. |15j . 
The brown line is the insertion work pWins produced by the FMT theory. Z = distance from 
the wall. Unit of length is taken to be the hard-sphere diameter a. The horizontal line is at 
p* — 0.5745 (a guide for the eye). The pseudo-density p pseu do(z) oscillates weakly around the 
bulk value and in opposite periodicity as compared to /?£ (z). The insertion work f3Wins also 
oscillates synchronously with p psC udo(z) but oppositely to p^\z). 



The wall potential w on the fluid particles (in the z-direction perpendicular to the wall) is 

fw(z) = oo, VzsCf, ^ 



I w(z) = 0, \fz > 



Thus, the closest approach of the center of a hard sphere to the wall is at z = a/2 (we shall use 
the hard sphere diameter a as unit of length below). For p* = 0.5745, the pure hard-sphere fluid 
excess chemical potential (3p,Q is 4.897, as calculated from the CS equation. The singlet density 
profile pw\z) has been obtained by Groot et al. [TS] via MC simulation and is displayed in fig ure 1 1 1 
(see the diamond symbols {>). The insertion work j3W{ ns {z) is identified as — Cw (z) (equation |L 



i.e. the (negative) 1-DCF, which is obtained from the MC data via the EL equation (equation ( 1.5 1. 
The insertion work is plotted as triangles (a). We remark that both the singlet density Pw (z) and 
the insertion work f3Wi ns (z) are "real" quantities, i.e. they are measurable quantities for the system 
of interest. We next formulate the pseudo-densities as mentioned earlier in section 1. We ask at 
what fictitious densities p pseu do(z) for a uniform HS fluid (which follows the CS equation), will we 
obtain a chemical potential value f3fi a that is equal to the MC nonuniform insertion work f3Wi ns (z) 

/?Mo = ^— ^ + Z' m = (3WUZ), (3.3) 

where n — |p psC udo (z) c 3 . In figure pi we show this mapping via a diagram. Note that Z' HS is 
the non-ideal compressibility factor from the CS equation. The pseudo-density p P scudo(z) thus 
calculated is also plotted (as asterisks *) in figure [l] We observe that the pseudo-density stays 
fairly flat and close to the bulk value of 0.5745, oscillating only weakly up and down with respect 
to the horizontal line. The oscillations are out of phase with respect (opposite in period) to the 
singlet density (z). The singlet density pw\z) has a contact value w 2.296 (z = cr/2). It 
oscillates vigorously with pronounced peaks and valleys. We note that the insertion work pW- ms {z) 
also oscillates with trends similar to the pseudo-density: when p pS cudo(z) is up, (3Wi ns (z) is also up, 
and when p pS eudo(<z) is down, /3Wi ns (z) is also down. 

A question naturally arises: how does this pseudo-density p psC udo(z) compare with the commonly 
used weighted densities, say, n a (z) in the fundamental measure theory (FMT) of Rosenfeld [27 . We 
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Figure 2. "Mapping" of properties between the non-uniform system and the uniform systems. 
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Figure 3. Comparison of the effective (pseudo) density p psou do (symbol *) with the weighted 
densities no (green line) and ns/V (blue line) from the FMT theory of Rosenfeld. Bulk hard 
sphere density p* — 0.5745. Symbol (^>) = the singlet density profile pw(z) from MC data of 
Groot et al. I15I. 



selected two weighted densities uq(z) (green line) and (z)/V (blue line) from FMT for comparison 
in figure [3] (there are in fact six weighted densities n a (z) in FMT (a — 0, 1, 2, 3, VI, V2) based 
on the geometries of hard spheres). The basis set of weights are composed of three functions: oj 2 
is based on the Dirac Delta function, thus generating n , u>3 is based on the Heaviside function, 
generating n 3 , and cuv2 is a vector function (see Rosenfeld [27J). The other linearly-dependent 
weighting functions are related to this basis set by scaling through spherical geometries. The first 
three weighted densities are related by the relations: uq(z) = ni(z)/R = 112(2) /S (R and S are the 
radius and surface area, respectively, of a hard sphere), while n 3 (z) is normalized by the spherical 
volume V : n 3 (z)/V. 

In figure [3] we notice that the magnitudes of both uq(z) and n^,(z)/V are close to the pseudo- 
density Ppseudo^) and far from the singlet density pw (z). There are visible differences between 
Ppsoudo (z) and no (z) as well as n-i(z)jV . This is not surprising, since these coarse-grained densities 
are derived from different theories (the FMT vs. the PDT). p psC udo(z) is close to no{z) only near the 
wall and at long ranges. The cusp in n (z) (due to the delta weighting) is missing from p psC udo( z ) ■ 
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Values of n 3 (z)/V are quite distinct from both Uq{z) and p P seudo(- z )- The observations also apply 
to other two weighted densities owing to scaling: n\{z)/R and ri2(z)/S. On the whole, the weighted 
densities and the pseudo-density are commensurate with each other (of similar magnitudes) even 
though they came from entirely different origins. 

To make further comparisons, we also calculate the insertion work arising from the FMT. The 
excess free energy F cx in FMT according to literature [27]) is expressed as: f3F cx = J dr$(n a (r}) 
where $ is a free energy density. This quantity is a function of the weighted densities n a and itself 
arises from a uniform fluid model (such as the compressibility equation from the Percus-Yevick 
equation |T7J (PY C ) or variations thereof [2S]). The FMT insertion work can be evaluated as the 
functional derivative [IE] SF CX /Spw\z), thus in Fourier space 



/3W ins{FMT) (k) = -^—uj a (k). 



(3.4) 



This FMT insertion work (W ins (FMT)) at p* = 0.5745 is plotted in figure [T] (the brown line). It 
quite closely follows the simulated Wi ns (MC) (triangles). This is to be expected, since both the 
FMT and the PDT-based theories are designed to match the nonuniform free energy. The FMT free 
energy lies somewhat above the PDT values, because the PY C solution for the chemical potential 
is higher than the MC value for hard spheres. 
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Figure 4. The effective (pseudo) density p pS eudo (symbol *; from PDT) and the insertion work 
PW-ms (symbol A; from MC) at bulk hard sphere density p* = 0.715 for the HS/HW system. 
The insertion work is scaled down 50% in order to fit in the graph. Symbol (<0) = the singlet 
density profile p^\z) from MC data of Groot et al. [15] . The brown line is the insertion work 
/3Wins produced by the FMT theory. The horizontal line is at p* = 0.715. The pseudo-density 
Ppseudo (z) oscillates weakly around the bulk value and in opposite periodicity as compared to 
Pw (z). The insertion work /3Wi ns also oscillates synchronously with p peC udo(z) but oppositely 
to p2\z). 



Similar calculations are made for HS/HW at a higher density p* — 0.715. Figure [4] shows 
the singlet density pw\z) (diamonds). The pseudo-density (asterisks) oscillates mildly. The FMT 
insertion work (brown line) is close to the MC insertion work j3W- lns (triangles) but overestimates 
it. The PY C values are higher than the MC values as already noted. (N.B. The values of pW- lns are 
multiplied by 0.5 (or 50%), in order to fit into the same graph). 

In figure |5] /9p SCU do (z) is compared with tiq(z) (green line) and n^(z)/V (blue line). Again, we 
see that at a few points, hq(z) and p pse udo(z) match closely, while n^(z)/V has larger differences 
from Pp S0U do (z) at this higher density. The comparisons show that the PDT approach is clearly 
different from the fundamental measure theory. FMT was constructed based on the PY C theory 
(or modifications thereof) in order to map the compressibility free energy (see also later improved 
versions [29 ) to the nonuniform free energy. In our approach, we use the PDT to map the CS 
free energy to the nonuniform free energy, and we use MC data to obtain the insertion work. The 
objectives are the same; but the theoretical bases are different. 
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Figure 5. Comparison of the effective (pseudo) 
densities no (green line) and ns/V (blue line) 
sphere density p* — 0.715. Symbol ((}) = the 
Groot et al. [Tg] . 



density p psC udo(z) (symbol *) with the weighted 
from the FMT theory of Rosenfeld. Bulk hard 
singlet density profile pw (z) from MC data of 



To further examine high-density hard spheres, we next examine the case p* = 0.813. The 
results are plotted in figure [6] Similar oscillatory behavior is observed for the pseudo-density and 
the insertion work (reduced to 20% in order to fit in the graph) . The phases of oscillations of the 
singlet pw\z) and p pS cudo(z) (and /3Wi ns as well) differ by a half period it. 




0.5 1 1.5 2 2.5 3 



Figure 6. The effective (pseudo) density p pscu do (symbol *; from PDT) and the insertion work 
PWvas (symbol A; from MC)) at bulk hard sphere density p* = 0.813 for the HS/HW system. 
The insertion work is scaled down to 20% in order to fit in the graph. Symbol (<0) = the singlet 
density profile Pw (z) from MC data of Groot et al. The horizontal line is at p* = 0.813. The 
pseudo-density p pS eudo oscillates weakly around the bulk value and in opposite periodicity as 
compared to p^\z). The insertion work f3W\n S also oscillates synchronously with p pS oudo(^) but 
oppositely to Pw\z). 



4. Construction of a closure theory 

It has been recognized that the DFT for the structures of nonuniform fluids should first be 
based on valid uniform-fluid properties as inputs. In particular, the structure of nonuniform fluids 
depends critically on the accuracy of the pair direct correlation functions 2-DCF of the uniform fluid 
that must be supplied beforehand. We have developed in the last decades an accurate theory |7] 
for the uniform hard sphere fluid: the zero-separation-theorem-based closure (ZSEP). We can thus 

(2) 

easily generate accurate 2-DCF Cq HS for use in the present study of nonuniform hard-sphere fluids. 
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We shall look at the five densities where the MC data are available for the HS/HW systems: 
namely p* = 0.5745, 0.715, 0.758, 0.813, and 0.9135. The ZSEP closure for the uniform bridge 
function -Bqhs is of the form [7] (caret indicates a function) 



Sohs(7*) = -^7* 2 



1-0 + 



«7 



(4.1) 



where a, 0, and £ are adjustable parameters, and 7* = 7 + p^. / is the Mayer factor of the repul- 
sive WCA (Weeks-Chandler- Andersen) potential. For details, see references in [7]. The adjustable 
parameters a, <j>, and £ arc determined using structural consistencies (the zero-separation theorem 
and the contact- value theorem) and thermodynamic consistencies (e.g., the pressure consistency, 
and Gibbs-Duhem relation). The uniform 2-DCF's C^ s obtained are depicted in figure [7J They 
have all satisfied closely the consistency relations imposed. For example, in the case of the high 
density p* — 0.90, the C^ s is compared with the MC-generated data [30J in figure [s] Diamond 




Figure 7. The pair direct correlation functions C^ s (r) of uniform hard spheres generated by 
the ZSEP closure [7 for five densities: p* = 0.5745, 0.715, 0, 758, 0.813, and 0.9135 (Curves from 
top to bottom, respectively). They all show a non-zero tail at r > 1 (a correct behavior for hard 
spheres) . 



CnHS 




Figure 8. The pair direct correlation functions CqSa(r) °f uniform hard spheres generated by 
the ZSEP closure [7| for the bulk density p* = 0.90. The lines are from ZSEP theory, and the 
symbols ({>) from MC data of Groot |30| . The inset shows the magnified view of the non-zero 



tail and discontinuity in C^%(r). At r = cr+, C ( ^ s (a + ) = 1.20 (from ZSEP), and 1.05 (from 



(2) 



MC). At r = a~, Cohs(°" _ ) = -3.965 (from ZSEP), and -4.12 (from MC). 



symbols denote the MC data, and the solid line is the 2-DCF from the ZSEP (4.1). We see very 



(2) 

close agreement of the two curves. The zero-value Cq HS of MC is —44.90 while the ZSEP value is 

(2) 

—47.70 (within 6%). In particular, the values CQ HS (r) derived from ZSEP for r > 1 are positive and 
non-zero. This is corroborated by the MC data. The positive peak C^ s (cr + ) = 1.2 from ZSEP, 
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while MC gives 1.05. This is in obvious contrast to the PY-based theories that would produce 
Cq HS (r) = for r > 1 (an erroneous result) . The same observation applies to all five present cases 
(figure [7]) . 

For the nonuniform HS/HW case, we developed a closure similar to ZSEP [3_: 



B w {"t w ) = -sgn( 7tu ) • Cll 



1-6' 



1 + a'j u 



(4.2) 



where a' , 6' , and £' are adjustable parameters for the nonuniform B w . sgn(.) is the sign function. 
j w is the indirect correlation function defined as 



7tt ,(f) = / dr'C ( V s (\r-r'\)5 Pw (r'). 



(4.3) 



We found from practice that "renormalization" of the indirect correlation is not needed for this 



nonuniform fluid. We also note that the expression (4.2 1 is modeled after the ZSEP equation (4.1 1 
for uniform fluids [7] . In earlier studies [2 [3] we have shown that this closure was effective for 
adsorption of high-density Lennard- Jones fluids. The similarity between the uniform closure and 
the nonuniform closure used represents an ongoing effort in the transfer of successful theories for 
uniform liquids to nonuniform liquids. 



Table 1. Parameters for the ZSEP closure (4.2 1 and wall sum rule for the HS/HW system 



Po 


0.5745 


0.715 


0.758 


0.813 


0.9135 


a' 


0.4 


0.4 


0.4 


0.5 


0.62 


4>' 


0.4 


0.4 


0.4 


0.4 


0.4 


C 


-0.38 


-0.28 


-0.22 


-0.20 


-0.13 




2.288 


4.319 


5.131 


6.598 


10.23 




2.29 


4.27 


5.15 


6.57 


10.28 



t Contact value (from ZSEP calculations). 
^Contact value (from hard-wall sum rule: p^}\a/2) 



P/kT) from Carnahan-Starling equation. 



Equations \1.6\ , (4.2 1, and (4.3 ) are coupled. The convolution was solved by standard numerical 
methods [3 in bipolar coordinates. Trapezoidal rule was used in integration. Picard's iterations with 
mixing (relaxation) on input and output iterates were employed. The grid of numerical integration 
was Ar = O.Oler, and total grid number N — 2048. Thus, the maximum z-distance reached is 
20.48cr. Cauchy's absolute convergence was enforced for the 7^,-function with convergence criterion 
S = 0.0001. There were three free parameters to be determined in equation (4.2): a', 6', and 
We needed three conditions. The first condition was the hard- wall sum rule: the contact density 
Pw\cr/2) = P/(kT). In practice, we set 0' to a default value of 0.4. a' and £' were adjusted in 
tandem until the wall rule was satisfied (see table [T]). In the future, other sum rules [3T] will be 
tested and incorporated. 

Figure [9] shows the result of application of ZSEP. The density varies from the moderate p* = 
0.5745 to the highest p* = 0.9135. The symbols are the MC data [15], the lines are from the 
ZSEP theory. We observe overall superb agreement for densities p* from 0.5745, 0.715, 0.758 to 
0.813. At p* = 0.813, our results compare well with previous accurate theoretical attempts (e.g. 
the White-Bear version of FMT [2S])- The ingredients for success are from (i) the high accuracy 
of the input uniform 2-DCF's; and (ii) the effectiveness of the ZSEP closure ( |4.2| . For the highest 
density p* = 0.9135 there is some deterioration of the theoretical curve at the second peak. The 
heights of the peaks (of ZSEP and MC) are similar, however the predicted location is displaced to 
a longer distance. A closer reading of the original MC source [15] reveals that at this high density 
(packing fraction r\ = 0.4783), there are clusters of molecules forming hexagonal structures for the 
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p. 



12 
10 

8 + : 

6 J> i> 



4 



0.5 



1.5 



0.9135 
0.81 30 
0.7580 
0.7150 
0.5745 



2.5 



Figure 9. The nonuniform singlet density functions pw\z) obtained from the new ZSEP-closure 
equation ( [4T2J for the HS/HW system at five densities p* = 0.5745, 0.715, 0.758, 0.813, and 
0.9135 (from bottom to top. Each curve shifted up by 0.5 units, respectively). Symbols = MC 
data from Groot et al. [15]. Lines = ZSEP results. We discern close agreement between the 
ZSEP curves and the MC curves except at the highest density (0.9135) where the second peak 
of the theoretical curve is displaced to a larger z- value. (See discussions in the text). 
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14 
12 



£ 10 
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0.5745 







0.5 



1.5 



2.5 



Figure 10. The insertion work generated by the new closure theory at p* = 0.5745, 0.813, and 
0.9135 (Curves from bottom to top, respectively). Lines = ZSEP theory; triangles = MC data. 
Close agreement is in evidence for p* — 0.5745 and 0.813. The deviations at p* = 0.9135 reflect 
the same cause: the two-dimensional solid-fluid transition of the surface layer of hard spheres at 
this high density. 



first layer of particles (figure [6] of |15j ) , and reduced normal-direction diffusion coefficient for the 
surface molecules: all indications of two-dimensional quasi-solid-fluid transition. Note that the fluid 
limit of homogenous hard spheres is at about 77 « 0.4948. For inhomogeneous hard spheres this 
may occur earlier. This puts stringent demand on a closure theory which is not intended to predict 
solid formation. 



Next, we calculate the insertion work via equation (2.141 using ZSEP-density pw\z) as input. 



The results are shown in figure |T0] for three densities p* — 0.5745, 0.8f3, and 0.9435. The agreement 
between the ZSEP (3W ins (lines) and MC /3W ins (triangles) is excellent for p* = 0.5745 and 0.843. 
For p* — 0.9435, similar discrepancies as for the singlet density pw\z) are in evidence. The 
displacement of the first minimum to a larger distance corresponds to the displacement of the 
second peak of pw\z). 
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5. Conclusions 



The two major questions that we set out to explore in this study are (i) how does the effective 
density based on the potential distribution theorem look like and how does it behave? (ii) Can 
we obtain the structure of hard spheres adsorbed on the hard wall through a closure theory for 
nonuniform fluids? We formulate two potential distribution theorems: one for the insertion work of 
the usual test particle that interacts with the wall, and the other for the test particle that ignores 
the wall. It is the latter insertion work that corresponds to the singlet direct correlation and as a 
consequence the second potential distribution theorem becomes equivalent to the Euler-Lagrange 
equation. 

The pseudo-densities p ps0 udo , due to their construction, are much smoother quantities than the 
singlet densities. They oscillate around the bulk value with periods out-of phase with respect to 
the nonuniform pwK 

The insertion work is derived from the MC data on HS/HW (Groot et al., 1987 [15]) via the 
Euler-Lagrange equation (with the identification f3W- ms = — Cw )■ Its behavior for several densities 
(p* = 0.5745, 0.715, and 0.813) is determined. j3W- lns oscillates in phase as the pseudo-density, but 
out of phase with respect to the singlet density pw ■ 

(2) 

The density functional theories need accurate uniform pair direct correlation functions Cq HS as 
input. Earlier we have formulated an accurate closure theory for the 2-DCF (the ZSEP theory) of 
uniform hard spheres. It is computationally cheap now to generate these functions: it took less than 
a minute of CPU time on a PC computer to produce an accurate DCF. Here we propose a new 
closure theory for the nonuniform hard spheres equation (4.2). Its form is inspired by the uniform 
fluid ZSEP equation. This bridge function is used in the star-function based density functional 
theory equation (1.6 1 to generate the singlet density profiles. These ZSEP densities are compared 
with the MC data (figure ^ at five densities p* = 0.5745, 0.715, 0,758, 0.813, and 0.9135. Except 
for the highest density, the ZSEP gives accurate nonuniform density profiles. This shows that 
the closure-based density functional theory can perform reasonably well for hard spheres when a 
suitable closure (bridge function) has been adopted. We shall have further developments on this in 
the future. 

These theoretical predictions for the singlet densities are used to calculate the insertion works 



(figure 10 1. They also compare favorably with the MC data. 



In summary, we have proposed and evaluated a new type of effective densities via the potential 
distribution theorem. We characterize their behavior for the hard spheres/hard wall system at 
moderate to high densities (up to p* — 0.9135). The free energies (the insertion works Wi ns ) of the 
nonuniform system are also calculated. They show oscillations opposite in phase to those of the 
singlet densities. We also construct a new closure-based density functional theory that can give 
accurate reproduction of the computer simulated densities and insertion works. 

The newly developed star-based density functional theory can be extended to fluids interacting 
with soft/attractive potentials. This task is under active investigation. 
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AflCop6i4m TBepflnx ccpep: CTpyKTypa Ta ecpeKTMBHa ryc~ri/iHa 
BiflnoBiflHO flo TeopeMM po3nofli/iy noTei-mia/iy 

n.n. /iR r. UemKa^ 

Biflfli.n xiMinHoT iH>KeHepiT Ta iH>KeHepiT MaTepiajiiB, YHiBepcMTeT KajiicpopHiT, 
rioMOHa, Ka/iicpopHia, CLUA 

LLlKcna cpi3HKH, y h iBepcnTeT iMeHi KBa3y.ny-HaTa.nb, CKOTCBmb, 
3209 niTepMapii46ypr, niBfleHHa AcppwKa 

Mm nponoHyeMO hobmm Tun ecpeKTU bh nx rycTMH, OTpwviaHi/ix 3 TeopeMW po3nofli/iy noTeHLn'a.ny. L|i rycTHHH 
Heo6xiflHi fljiiR Biflo6pa>KeHH!R Bi/ibHoT eHepriT oflHopiflHoY piflMHi/i y BijibHy eHepriio HeoflHopiflHoY piflMHi/i. Teo- 
peMa po3noflijiy noTeHLiia/iy flae po6oTy, Heo6xiflHy p,nn ycTpoMJinHHa npo6HoT HacTn h km y cucTeiuiy Mo.neKyn, 
Ha siKy flie 30BHiujHi£i noTeHL4ia/i. L(a po6oTa VK; ns Mox<e 6yTH OTpi/iMa Ha 3 CMMynnLn'T MoHTe-Kap.no (MK) 
(HanpwKJiafl, niflxoflOM TecTOBoT HacTi/iHKW BiflOMa) a6o 3 aHa/iiTWHHo'f TeopiT. riceBflorycTWHH e no6yflOBaHi 
TaK, mo KO.ni/i "fx 3HaneHHn niflCTaB/iaKDTbca b piBH^HHH CTaHy oflHopifiHoY piflHHH (HanpuKJiafl, piBHflHHfl fl/ifl 
xiMinHMX noTeHLiiajiiB KapHaraHa-CTap^iHra CHCTeMU TBepflwx ccpep), to BiflTBoptoeTbca po6oTa ycTpoMJiflH- 
hh HacTWHKH b HeoflHopiflHy pifli/iHy, OTpi/iMa Ha 3 ci/iMynsmiT MK. Mh floc/iifl>KyeMO noBefliHKy nceBflorycTHHH 
p/in CMCTeMM "TBepfli ccpepi/i"-"TBepfla CTiHKa" npn cepeflHix Ta bwcokmx rycTHHax (Bifl p* = 0.5745 flo 
0.9135). fl/ia Liboro bh Kopi/icTOByKDTbcs pe3y/ibTaTM MoHTe-Kapjio TpooTa Ta cniBaBTopiB. riceBflorycTHHH 
fleMOHCTpyKDTb ocL4MiiK>KDHy noBefliHKy 3 npoTWJie>KHOK) cp330KD flo oflHonacTWH kobi/ix rycTWH. Mh t3ko>k npo- 
noHyeMO HOBy Teopiro cpyHKL4ioHa.ny rycTi/iHH Ha ocHOBi 3aMMKaHHfl (reopiio dpyHKujonany rycTMHi/i 3ipKOBoY 
cpyHKLiiY), U40 MO>Ke tohho o n 1/1 cy Ba t 1/1 npocpmi rycTHHH i po6oTy ycTpoM/iaHHa. ToHHicTb TeopiT nepeBipaeTbOR 
fl/i!R fleKijibKox BwnaflKiB aflcop6i4iT TBepflwx ccpep. 

K/iK>HOBi c/iosa: reope/wa po3nop,)ny noreHLtiajiy, TBepfli cdpepiA, ap,cop6u,\n, ecpeKTHBHa rycTWHa, 
CTpyKTypa, 'mTerpanbHe piBHHHHn, 33MtiKaHH& 
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